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We present lattice calculations in QCD using a variant of Kaplan 

fermions [l],^] which retain the continuum SU(N)l x SU(N)r chiral symmetry 

j-C". on the lattice in the limit of an infinite extra dimension. In particular, we 
ON. 

, show that the pion mass and the four quark matrix element related to K q — Kq 

| mixing have the expected behavior in the chiral limit, even on lattices with 

o _ 

^ | I modest extent in the extra dimension, e. g. N s = 10. 

m 

>. 

A fundamental property of the QCD Lagrangian is chiral symmetry, or invariance under 
independent global flavor rotations of the left handed and right handed quark fields. Through 
■ the familiar Goldstone theorem, the spontaneous breakdown of this symmetry generates 



massless excitations, and hence chiral symmetry is important to the dynamics of QCD. In 
reality the symmetry is also explicitly broken by small quark masses, so the excitations 
are not massless, but very light compared to masses of other hadrons. Until now, lattice 
simulations of QCD have only accounted for this symmetry approximately or not at all, 
relying on the continuum limit to restore it exactly. Since most simulations are far from the 
continuum limit, this is not satisfactory. 

The two popular discretizations of the Dirac action used to date for lattice simulations 
of QCD are the Wilson quark action and the Kogut-Susskind action. For non-zero lattice 
spacing, each breaks the continuum SU(N)l x SU(N)r chiral symmetry. The Wilson action 
includes a mass-like term, and the Lagrangian quark mass renormalizes additively. In this 
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case, the chiral symmetry breaking is severe, and fine tuning of the bare parameters is 
required to obtain a massless pion, which defines the chiral limit. Indeed, these Goldstone 
bosons may not even be related to chiral symmetry in the nonperturbative regime, but 
instead to the spontaneous breakdown of parity-flavor symmetry ||. In addition, Wilson 
fermions lead to mixing of operators of differing chirality which is absent in the continuum. 
Physical results depend crucially on removal of these effects. This has been a long-standing 
problem in lattice calculations of weak matrix elements with light hadrons f3]|| that use 
Wilson quarks. Weak coupling perturbation theory is unable to accurately match the lattice 
operators to their continuum counterparts. Various methods have been developed over the 
years to deal with this problem, and so far none is very satisfactory [§-|9|]. 

Kogut-Susskind fermions retain an exact U(l) x U(l) remnant of the continuum 



SU(N)l x SU(N)r chiral symmetry QTO | , so the bare quark mass renormalizes multiplica- 
tively, and the chiral limit is obtained as in the continuum, by taking the quark mass to 
zero. Because of this remnant chiral symmetry, Kogut-Susskind quarks are preferred for 
light hadron matrix element calculations and simulations of QCD at non-zero temperature 
where it is now widely believed that QCD undergoes a phase transition that restores chiral 
symmetry. This phase transition happens at a critical temperature of roughly 150 MeV, 
or approximately the mass of the pion. Thus, the correct number of pions is important 
to describe the dynamics of the phase transition. However, at finite lattice spacing the 
Kogut-Susskind action has only one Goldstone pion instead of the N 2 — 1 in the continuum 
theory. 

For these reasons it is desirable to use a discretization of the QCD action that preserves 
chiral symmetry. Some time ago Kaplan [|IJ proposed a discretization of the Dirac action 
for chiral fermions that avoids the notorious fermion doubling problem and hence could be 
used to formulate a nonperturbative chiral gauge theory. The impossibility of formulating 
such a theory that is local and gauge invariant is summarized by the Nielsen-Ninomiya no-go 
theorem [ IT| . For each chiral fermion that one starts with in the continuum, a doubler of the 



opposite handedness appears that renders the lattice theory vector-like. This phenomenon is 



quite general and occurs already in free field theory. However, for free Wilson-like fermions 
in d+1 dimensions, Kaplan showed that including a mass term with a defect along the extra 
dimension results in a single chiral fermion which is bound to the d-dimensional defect. Thus 
a low energy effective theory of chiral fermions may be described on the lower dimensional 
subspace of the defect. On a finite lattice with periodic boundary conditions, an anti-defect 
appears with a bound chiral fermion of opposite chirality, so again the theory is vector-like. 
However, if the extent of the extra dimension is large, the fermions have exponentially small 
overlap and do not mix; one is free to study the lattice chiral fermion on either wall. Sadly, 
when the chiral symmetry is gauged, the fermions mix, or if only one is gauged, new particles 



appear in the spectrum, and the chiral properties are probably lost fl2fl . But, see also the 



overlap formalism of Ref. [13 



The above considerations are irrelevant for QCD which is a vector gauge theory, i. e. right 
and left handed fermions couple to the gauge fields with equal strength. The important point 
is that Kaplan's fermions retain the full SU(N)l x SU(N)r continuum chiral symmetry on 
the lattice which makes them an attractive alternative for simulating QCD. Actually, at tree 
level one can show that the symmetry is broken explicitly by terms which are exponentially 
small in the size of the extra dimension [Q. The main purpose of this work is to verify 
that these terms remain small in nonperturbative simulations and do not spoil the chiral 
behavior of observables. 

Below we describe calculations in QCD with a variant of Kaplan's original fermions that 
has been proposed by Shamir p,l4|. In this version only half of the lattice in the extra 



dimension is used. The chiral fermions, one left handed and one right handed, reside on 
opposite boundaries of the five dimensional space, and they interact with four dimensional 
gauge field configurations which are the same on each slice in the five dimensional world. The 
gauge links in the fifth dimension are set to the identity except for the links which connect 
the boundaries. Those links explicitly couple the two chiral fermions with a strength equal 
to — m. In perturbation theory it was shown in Ref. that m is proportional to the bare 
quark mass of the vector fermion constructed from the left handed and right handed quarks 



on the opposite walls of the fifth dimension, m q = mM(2 — M). Here M is the ordinary 
Dirac mass in lattice units for a five dimensional fermion. Then in Ref [I5[] it was shown 
that as the number of lattice sites in the extra dimension, N s , is sent to infinity, the axial 
currents of domain wall fermions satisfy nonperturbative Ward identities. Thus the chiral 
limit is obtained by simply setting m = 0, unlike in the case of regular Wilson fermions. 

For details on Kaplan fermions and the variant of Shamir we refer the reader to 
Refs. 0,|,15]; here we simply state the quark action for one of the flavors [[15 . 
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where ip is a five dimensional Dirac fermion. The U x li are the usual four dimensional gauge 
links and are elements of the color group SU(3). Note the Wilson term which removes the 
doublers is added to the action with minus the conventional sign, x and y denote four 
dimensional coordinates while s and s' denote coordinates along the extra fifth dimension. 
A^ s is the number of sites in this dimension. The simplest choice for operators which create 
and destroy light four dimensional quarks are |I5[ 



(I+75) (l-7s) , 
Qx = Wxfi h y>x,N.-i, 

_ j (1 + 75) , T (1 - 75) , , 

In fact, any operators localized near the boundaries that have finite overlap with the bound- 
ary states will do as interpolating operators for the quark states. 

As a first step in implementing domain wall quarks for QCD, we measure observables 
constructed from the above operators on existing configurations. The gauge field configura- 
tions that we have used are a set of 20 lattices generated with two flavors of Kogut-Susskind 



dynamical quarks at gauge coupling 6/g 2 = 5.7 and bare quark mass rnxs — 0.01 [fL6| . The 
lattice size for these configurations is 16 3 x 32. In all of our calculations, the Dirac mass M 
was set to 1.7. 

For free quarks one must take < M < 1 to obtain massless states with strictly positive 
Green's functions [Q] . For 1 < M < 2 there exists a massless mode, but its Green's function 
oscillates like (1 — M) s . This free field condition is renormalized at strong coupling as we 
found only heavy states for M < 1. This can be seen from perturbative arguments || and a 
simple mean field treatment which gives 4 — 4Uq < M < 6 — AUq ||17|| . Here Uq is the fourth 
root of the plaquette expectation value, or the mean field value of the gauge link. The mean 
field propagator then falls off like (5 — (M + 4Uq)) s . On the above configurations Uq ~ 0.577, 
yielding M ~ 1.52 for the most rapid decay of the propagator in the extra dimension, and 
thus the smallest overlap of the boundary states. The above mean field result agrees fairly 
well with the "optimal" value of M — 1.7 which we determined by observing the exponential 
decay of the quark propagator in the extra dimension. As M departs from this optimal value, 
the extra dimension must be increased to avoid possible overlap of the massless modes. While 
massless pions may still exist for M < 1 at strong coupling, we did not observe them due 
to the relatively small extent of the lattice in the extra dimension. Thus, for efficient lattice 
simulations with domain wall quarks it is important to determine the optimal value of M. 
In this study we set the number of sites in the fifth direction to N s = 4 and 10 to study the 
effects of finite size in the extra dimension on the chiral symmetry properties of observables. 
Below we will see that for M = 1.7, N s = 10 is sufficiently large to suppress the overlap of 
the light modes. 

The simplest test of chiral symmetry (breaking) for domain wall fermions is to measure 
the pseudo-scalar two-point correlation function and extract the mass. Chiral perturbation 
theory requires this mass to go to zero like \/m, or m% cx m. In Fig. [I] we show m% versus m 
for the two values of N s . The masses were extracted from a single particle covariant fit to 
the pseudo-scalar correlator at large Euclidean time separation from the source (see below). 
At N s = 10, the lightest two points linearly extrapolate to 0.0002 ± 0.0160 while for N s = 4 



the results clearly miss the origin to the left. It appears that N s = 10 is large enough to 
suppress the explicit chiral symmetry breaking effects due to the overlap of the two chiral 
modes. The N s = 4 case shows that the approach to the N s — > oo limit is from above which 
is expected since the Wilson term is added to the action with minus the conventional sign. 

For a more stringent test of the chiral symmetry properties of domain wall quarks we turn 
to the calculation of hadronic matrix elements, which as mentioned above is problematic with 
conventional lattice operators. The canonical example is the Kq — Kq mixing matrix element, 
Mn, required to determine the CP violating phase of the CKM matrix. Explicitly, Mil = 
(K\[s , ~f u (l — '-)5)d\ 2 \K) and is expected to vanish as m 2 K in the chiral limit. However, as we 
have mentioned, it is well known that the matrix element of this operator constructed from 
Wilson quarks does not vanish in the chiral limit. Only by fine tuning a linear combination of 
all four quark operators with differing Dirac structures does one obtain the correct behavior. 
To date, no satisfactory method of fine tuning has been found to solve this problem 0-0, 
so lattice operators which avoid this mixing are much preferred, cf. Kogut-Susskind quarks 
and now domain wall quarks. 

In Fig. we show (^|[s7„(l — 7 5 )c/] 2 |K)/|(0|s7 5 (i|_ft')| 2 versus m for both values of N s 
and the usual Wilson quark operator whose mixing coefficients were determined using lattice 



weak coupling perturbation theory |18|J19(| with a boosted coupling |20j. This ratio has the 
same mass dependence as Mll in the chiral limit since the denominator is just proportional 
to the square of the kaon decay constant. From Fig. ^| it is clear that the N s = 10 result 
exhibits the required behavior in the chiral limit: the values at the two smallest quark 
masses extrapolate linearly to .002 ± .015 at m = 0. On the other hand, the N s = 4 domain 
wall quark and Wilson quark results significantly miss the origin. We note that while the 
errors are considerably larger for the Wilson case, the results at each quark mass are highly 
correlated and a fit which accounts for this does not go near the origin. Indeed, previous 
lattice calculations have shown that increased statistics do not remedy the problem. Also 
note that the domain wall results again approach the large N s limit from above. 

The ratios in Fig. |2] corresponding to domain wall quarks are determined with the follow- 
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ing procedure. First we calculate five dimensional propagators from wall sources at times 
and 31 (the ends of the lattice) on each boundary in the fifth dimension. The quark prop- 
agators are then constructed using Eq. (0). Because of the chirality projection operators in 
Eq. (0), only two source spins on each boundary are required. We find that the added cost 
of inverting the five dimensional fermion matrix scales linearly with N s . Next, these propa- 
gators were contracted using point sinks to form three correlation functions on each gauge 
field configuration, a two-point pseudo-scalar correlation function from each wall source and 
a three-point correlation function from the same two wall sources. Each correlation function 
is averaged over configurations. The ratio is found by dividing the three-point correlation 
function by the product of the pseudo-scalar correlators, all three being evaluated at the 
same sink time slice. For large time separations between the sink time slice and the two 
wall sources, the result is a constant. Our results were averaged over sink time slices 10 
to 20 where this constant plateau was evident. The errors displayed in Fig. (Q), which are 
statistical, are calculated from a jackknife procedure. In fact, simultaneous covariant fits to 
all three correlation functions, which have good x 2 , give similar results. The Wilson quark 
results were calculated in similar fashion, but with point sources. 

A similar ratio as shown in Fig. |2], but with (0\sjsd\K) replaced by (0|s757 d|i^), gives 
the parameter Bk which chiral perturbation theory predicts to be non-zero in the chiral limit. 
We are in the process of calculating this important phenomenological quantity. Preliminary 
measurements indicate that the ratio depends weakly on m and is non-zero as m — > 0. The 
numerical value of Bk is scale dependent (as is Mll) and thus requires a renormalization 
factor. In addition, finite lattice spacing effects must be studied by simulating at several 
couplings to understand the continuum limit of Bk- Thus much work needs to be done 
before a meaningful value can be given. 

The above results show that domain wall quarks exhibit excellent chiral behavior, even 
at g 2 ~ 1. In particular the bare quark mass appears to receive only multiplicative renormal- 
ization, i.e. no fine tuning of the couplings is required to reach the chiral limit, in contrast 
to the case with Wilson quarks. Strictly speaking this is true only in the limit N s — > oo; 
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however, residual chiral symmetry breaking effects coming from finite N s , which were not 
detectable in our calculations with N s = 10, are in principle removable by simply increasing 
N s . These encouraging results also show that domain wall quarks are a vast improvement 
over Wilson quarks for matrix element calculations. Again, there is no fine tuning or cum- 
bersome operator subtractions necessary to obtain the correct chiral behavior. Similarly, 
they are likely to prove useful for other QCD calculations where chiral symmetry is crucial, 
e. g. the QCD phase transition for non-zero temperature where current simulations suffer 
from explicit chiral symmetry breaking effects. Finally, we re-emphasize that domain wall 
quarks at N s = 10 retain the full chiral symmetry of continuum QCD to a remarkable degree, 
unlike either Wilson or Kogut-Susskind quarks. 

We thank Y. Shamir and M. Creutz for suggesting the use of domain wall fermions 
for matrix element calculations and M. Creutz for helpful discussions. Our domain wall 



fermion code relies heavily on the four dimensional MILC code [^TJ, which we gratefully 
acknowledge. This research was supported by US DOE grant DE-AC0276CH0016. The 
numerical computations were carried out on the NERSC T3D and the SDSC Paragon. 
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FIG. 1. The pion mass squared as function of m. m is proportional to the quark mass. The last 
two points for N s = 10 (octagons) extrapolate linearly to zero well within statistical errors. The 
squares denote results for N s = 4. 
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FIG. 2. The ratio of the four quark matrix element for Kq — Kq mixing to the square of the 
pseudo-scalar density matrix element, calculated with domain wall fermions (octagons(iV s = 10) 
and squares(iV s = 4)). The N s = 10 curve exhibits the correct behavior in the chiral limit. Also 
shown is the result using the same gauge field configurations for Wilson quarks (crosses) which 
extrapolates to zero far from m = (note that for Wilson quarks the quark mass is defined as 
the difference of the inverse quark hopping parameter with the inverse critical hopping parameter, 
m = i^" 1 - K" 1 )). 
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